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An absolute lower bound on the number of templates needed to keep the fitting factor above a 
prescribed minimal value F in correlator-bank detection of (newtonian) gravitational wave chirps 
from unknown inspiraling compact binary stars is derived, resorting to the theory of quasi band- 
limited functions in the L°° norm. An explicit nearly-minimum redundant cardinal-interpolation 
formula for the (reduced, noncoherent) correlator is introduced. Its computational burden and 
statistical properties are compared to those of the plain lattice of (reduced, noncoherent) correlators, 
for the same V. Extension to post-newtonian models is outlined. 
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INTRODUCTION 

The detection of gravitational wave (henceforth GW) chirps from unknown inspiraling compact binary sources 
(henceforth CBS) is a primary goal for the early operation of broadband interferometric detectors, including TAMA300 
§, GEO600 @, the two LIGOs f| and VIRGO (|, in view of the sizeable expected rate of observable events @. 

For additive gaussian stationary noise, the correlator-bank threshold-detector is the optimal one, yielding the 
smallest false-dismissal probability, at any fixed false-alarm probability and signal to noise ratio Q . 

The issues of optimum template parametrization and placement, and the related computational burden have been 
discussed by several Authors JtJ- |l3| . A lucid account of the main relevant landmarks is given in 

Curiously, the question of possible efficient interpolation among the correlators has been left yet unsolved (li) . 

In this paper we set up and test an efficient interpolated representation of the (reduced, noncoherent) correlator for 
the simplest paradigm case of newtonian chirps. The proposed representation JT^ j is proven to get close to the absolute 
minimum template density required by a prescribed minimal-match condition, which follows from the theory of quasi- 
bandlimited (henceforth q-BL) functions in the L°° norm. The statistical performance of the proposed representation 
are shown to be essentially equivalent to those of the plain lattice. 

This paper is accordingly organized as follows. In Section I we recall a number of relevant concepts and results. 
In Section II we review the design of the plain template-bank for the simplest (newtonian) case, and discuss its 
statistical detection properties. In Section III we briefly introduce q-BL functions and cardinal expansions, and 
derive the proposed representations. In Section IV we compare the computational burden and the statistical detec- 
tion/estimation properties of the cardinal-interpolated (reduced, noncoherent) correlator lattice to those of the plain 
lattice. Conclusions follow under Section V. 



I. BACKGROUND 

In this section we resume a number of well known concepts relevant to CBS chirp detection, and introduce the 
notation. 



A. Noncoherent Correlator. Deflection and SNR 

Detecting GW chirps from unknown inspiraling CBS requires the computation of a suitable set of non-coherent 
correlators (henceforth NCC) Jig : 
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where (fmf, fsup) is the useful antenna spectral window, A(f) = S(f) + N(f) are the noise corrupted (spectral) data, 
resulting from the superposition of a (possibly null) signal S(f) and a realization N(f) of the antenna noise (assumed 
gaussian and stationary), T(/) is an element of a suitable set of unit-norm chirp- templates such that 
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and II(/) is the (one-sided) antenna noise power spectral density (henceforth PSD). The random variables c have 
Ricean probability densities B, 



w(c) = cexp I - — I I Q (cd) , 

where io(") is the modified Bessel function of first kind and zero order, and 
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is the deflection obtained using T. The moments of (1.3) can be written in terms of Kummer's confluent hypergeometric 
function [B5i: 



For d ~ 5 the riceans (1.3) merge into gaussians, with: 

S[c]«d, var[c]wl. 

The deflection attains its maximum value iff the template T is matched to the signal, viz. 

S(f) 



T(f) = 



\\S(f)\ 



yielding: 
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where SNR is the intrinsic signal-to-noise ratio 
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B. Chirp Templates 

The stationary phase principle (see JlTj for a thorough discussion of its validity) can be used to show that the 
asymptotic principal part (T^| of a general [^9|, |20|, reduced |2l[] post newtonian (henceforth PN) chirp can be 
written dill: 



S(f; 4> c , T c , g) = Af- V& exp [j [2nfT c - C + V(/, 6] } , 



(1.9) 



where A is a constant (real, unknown) amplitude factor, T c is the (fiducial) coalescency time [ p2[ , </> c is the template 
phase at t = T C , and ^ represents the remaining intrinsic source parameters [ p6| . 

Equation (1.9) is used to construct the needed chirp templates. A further suffix T will be used to label the template 
parameters At, 4> Ct , T Ct , and £t- 



All template amplitudes At will be chosen so as to comply with the normalization condition (1.2), viz.: 
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C. Maximum Likelihood Criterion. Fitting Factor 



Equations fll.q ) imply that (under the assumption of a uniform distribution of the unknown source parameters) 
the largest correlator will most likely correspond to the special template yielding the largest deflection (maximum 
likelihood (ML) estimation criterion ||). 

Data analysis for detecting chirps reduces thus to the following. Given the (spectral) noisy data, and a set 
(lattice) of templates, suitably covering the chirp parameter space, the corresponding (noncoherent) correlators 
{cfejfc = 1,2,..., N} are computed. The largest among these correlators is used as a detection statistic viz., 
whenever this latter exceeds a suitable threshold, set by the prescribed false-alarm probability (surveillance strategy 
H ) , a signal is declared to have been observed p4| ] , and the corresponding template is taken as the most likely estimate 
of the observed signal [||. 

It is convenient to measure the goodness of fit between a given signal S(f) and the best available template in the set 
in terms of the so called fitting factor p7| : 



FF = max . 

k SNR 



(1.11) 



The set of templates should be constructed in such a way that for any admissible signal, 



ff > r, 



(1.12) 



where 1 — T 3 gauges the fraction of potentially observable sources which could be lost as an effect of template mismatch 



D. The Reduced Correlator 



In view of eq. (l.S), the noncoherent correlator (1.1) can be written: 
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Equation ( 1.13 ) is formally (but for inessential factors) the absolute value of the (/ — » T CT ) Fourier transform of the 
(complex-valued) function: 



A(/)T*(/; 0,0,0 , , . . , 

K(f) = <( n(j) ' •"«/ — j - j su p 

0, / < /in/) / > fsup- 
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Maximizing the noncoherent correlator (1.13) w.r.t. T CT1 is thus equivalent to taking the largest absolute value of the 
(/ — » T CT ) Fourier transform of (1.14). The resulting reduced correlator Q will be denoted with a capital letter, viz.: 



C = sup c 

TV. 



(1.15) 



E. The Newtonian Deflection and the Match 



For illustrative purposes, in this paper we shall restrict to the simplest Newtonian (OPN) signals and templates. 
The OPN function Vt(/;C) in (|T|) reads: 

^(/) = I |(^) _5/3 ^ 5/3 /- 5/3 ) (1.16) 
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where M.t is the template chirp-mass. It is convenient to introduce the following dimensionless variables and param- 
eters: 

f - M 

f = T-> = f inf (T c -T CT ), M = —, 

Jinf Mq 



-5/3 



128 V c° 

where Mq is the solar mass, so as to recast the deflection (|1.4[) into the form: 



(1.17) 
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It is also useful to introduce the reduced deflection: 
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D(A) = maxd(A,8), 



and the related normalized |3ll functions: 
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The function -D(A) is known as the (newtonian) match. The quantity T in ( 1.12 ) is accordingly also named the 
minimal match [ 54[ . 

The functions J and Z) are displayed in Fig.s 1 and 2, respectively, for the special case of a LIGO-like noise PSD, 
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with /o = 300-ffz (n is a constant of no concern to us here), and for a spectral window with: 

finf = 40Hz, f sup = 400Hz. (1.22) 
The value Q ma x of 8 which maximizes d(A, 8) in a neighbourhood of A = is shown in Fig. 3 as a function of A. 



II. THE PLAIN (NEWTONIAN) LATTICE 

Given a range [M m i n , M max ] of allowed source chirp masses, let the set of template chirp-masses be: 
MT 5/3 = M-% 3 , M- 5 + { 3 = M-^ 3 + 6 L , n=l,2,...,N L - 1, 
where Sl is the lattice-spacing, and 



N L = 



is the total number of templates. Obviously A can take only the discrete values: 

A„ =M; 5/3 ~M- 5/3 , n = l,2,...,N L . 



(2.1) 



(2.2) 



(2.3) 
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A. Lattice Design 



Given A4 S , the fitting factor of the lattice is: 

FF = maxL>(A„) = D (min |A„| ) . (2.4) 



The minimal-match condition (1.12) should be obviously enforced in the worst case, where: 



yielding: 



M-^ = M- q b l*+ 5 -±, (2.5) 



Equation (2.6) uniquely determines the lattice spacing 8l, and hence via (2.2) also the total number of templates 



B. False Alarm and False Dismissal Probabilities 

The statistical distribution of the lattice detection statistic 

maxCfc = maxcf, (2.7) 

k h 

is not known in exact form, and one should resort to numerical simulations aided by intuition to compute the lattice 
false- alarm and false-dismissal probabilities. 

The detection threshold 7 is determined from the prescribed (tolerated) false alarm probability a, by solving the 
equation: 

a = prob[3 k:C k > -y\SNR = 0] = 1 - prob[V h, c h < j\SNR = 0]. (2.8) 



The joint probability in ( |2.8| ) is difficult to compute, since the Ch are not statistically independent in general. For 
most practical purposes, a decent approximation is: 

prob[V h, c h < j\SNR = 0] = (1 - e~"<'' '' /2 ) M , (2.9) 

which would be appropriate if the Ch were a collection of M independent ricean random variables. Numerical experi- 
ments suggest an almost linear dependence of M on the total number of NCC used [63] . 
The probability of false dismissal of a signal with SNR ^ is: 

/3(7,S , 7Vi?)=prob[maxC* fc <7|S , A^i?^0] = prob[Wi, Cfl < 7|SWi?^0[. (2.10) 

k 



A simple (conservative) approximation of ( 2.10) ) is |l^| , f{3| 



(3(7, SNR) « prob[C_<7,C+<7|5iVi? ^ 0], (2.11) 
where C_, denote the reduced correlators corresponding to the nearest-neighbouring templates, with chirp masses 



M.± such that Ai s G [A4-, M.+]. Under the same (reasonably large SNR) assumptions leading to (2.11), the involved 
joint probability density can be approximated by a gaussian bivariate [p4|: 



exp 



(C--d-) 2 + (C+-d+) 2 -2R(C--d^){C + -d + ) 
2(1 — R 2 ) 



(2.12) 



2 7 r(l-i? 2 ) 1 /2 
where: 

d± = d[A±,e max (A±)], R^d[S L ,e max (A + )-Q max (A-)}, (2.13) 
and the function Q m ax(') has been defined in Sect. II. F and shown in Fig. 3. 
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Letting: 



MJ 5/3 = M- 5/3 + f]S L , fjG[0,l[ 



(2.14) 



the false dismissal probability ( [2.11 ) is obviously a function of ?/. Within the limits of validity of (2.8) to ( 2.1 3| ) , for a 
fixed SNR and a prescribed a, the following qualitative dependence of (3 on the lattice spacing 8l is observed. In a 
neighbourhood of n = 0.5, the false dismissal probability is reduced by reducing the spacing St among the templates. 
On the other hand, in a neighbourhood of r\ = 0, reducing the template spacing produces an increase of /3. This is 
due to the dominant effect of the parallel increase of 7, needed to keep a unchanged. 



A judicious tradeoff should be obviously sought, to choice a value of T, and hence of 5l, via (2.6), which minimizes, 
e.g., the average value of fi w.r.t 77, under the assumption of a uniform distribution of the source chirp mass. 



III. A NEARLY MINIMUM REDUNDANT INTERPOLATED LATTICE 

This section contains the main new results. A short introduction to the theory of q-BL functions is included to 
make the paper self-contained. 

A. q-BL Functions and Cardinal Expansions 



A function / : x € [a, b] — ► 1Z is q-BL in the L°° norm iff 



3 7 ,B c eU+ : sup \f(x)-f B (x)\ = exp[-i(B - B e )], 

x£[a,b] 



where: 



W(x) 



1, M < 1 

0, \x\ > 1. 



(3.1) 



/b(x) being obtained by taking the inverse Fourier transform of the spectrum of f(x) chopped at |y| > B, viz.: 

f B (x) = F~i x [W (|) ■ T x ^ v [f(x)}} : (3.2) 



(3.3) 



For a strictly bandlimited function f(x), whose spectrum vanishes identically outside [— B, B], eq. (3.2) provides an 
exact interpolating representation known as cardinal expansion p& , p7l : 



fsix) = ^ fi x n) sine ~(x-x n ) 



x n+1 -x„ =6, S = 



(3.4) 



where sinc(a;) = sin(x)/x. For a q-BL function on the other hand, one can prove that eq. (3.4), while reproducing 
exactly f{x) at x = x^, k 6 N, satisfies eq. (3.1), i.e. that 



Ve > 0, 3 B 



sup \f(x) - f B (x)\ < e. 

.f'G [ff .fr] 



(3.5) 



Equation (3.4) is an approximate sample-interpolating representation, where the sample density S^ 1 = 2B depends 
on the prescribed L°° appr oxim ation error e. It is important to note that the exponential decay of the error in (3.1) 
implies that reducing e in (3.5) by orders of magnitude does not change the order of magnitude of B. 



Usually one needs to compute f(x) in a finite interval [a, b] including only 

n — '^ b - a y 



(3.6) 



samples 38|. However, using eq. (3.4) to compute f{x) in [a, b] requires, in principle, knowledge of an infinite number 
of samples outside the interval of interest. This limitation can be circumvented by using generalized (economized) 
cardinal expansions. These expansions have the general form pi: 
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f( X )= X! /(^n) sinC Tj( x ~ x n) 6(X — X n ), 



x n+1 -x n =S', S' = (2 x Br\ X >1 (3-7) 

where 0(x) is a suitable windowing function such that: 

0(0) = 1, ^ [0(s)] = 0, V £ > (* - X >1. (3.8) 

The expansion (|3.7|) is nothing but the std. cardinal expansion of the function f(x)9(u — x), whose bandwidth under 
the assumptions ( ft.q ) is B 1 — \B, viz.: 

oo 

f(x)d(u - x) = ^ f{ x n)0(u - x n ) sine — (x - x n ) 

n— — oo 

1 

x n +i - x n =6', 5' 



(3.9) 



evaluated at u = x. The Fourier spectrum of f{x)0{u — x) is a smoothed version of the plain spectrum of f(x)] a 
judicious choice of 6{x) can thus make in principle, the decay rate of f(x n )0(u — x n ) as \u — x n \ — > oo as fast as desired 

The Knab window function WOl : 



sinh < nP(l — x X ) 



6(x) = ifp(i) := 



1- 



\P5> 
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(3.10) 



sinh[ 7 rP(l-x- 1 )] 



satisfies all constraints ( |3.8| ) and is essentially confined [|4l| in |x| < PS'. This allows to truncate (3/7) at a; — a;„| « P<5', 
so that for any given x, only « 2P sam ples symme trical ly placed around x are essentially needed to reconstruct f{x). 
The error resulting from truncation of (3/7) with ( 3.1C ) at \x — x n \ — P5' has been discussed in |^o|. A simple (and 
conservative) upper bound is given by: 



^ f(x n ) sine — (a;-x„) /Yp(a 

x-x„|>P5' 



< 



< 



M . . 

smh[7rP(l -x x )] x6[ n ,i,] 



Usually one enforces the condition: 



M 



e' < e, 



sinh[7rP(l - x -1 )] 

where e is the prescribed L°° error in (U). Equation (|t|) can be solved to express P as a function of x> 

sinh^ 1 (M/e') 



^(i-x" 1 ) 



The total number of samples 



N T = 



(b-a) 



2P = 



(3.H) 



(3-12) 



(3.13) 



(3.14) 



needed to represent f(x) in [a, b] using ( |3.9| ) and ( 3.10 ), within a prescribed L°° error and under the constraint ( [3.12 ), 
can be accordingly minimized by letting: 



X = l 



28 sinh _1 (M/e') 
Tr(b — a) 



1/2 



(3.15) 
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B. The Cardinal-Interpolated Newtonian Match 



The match D(A) is a q-BL function in the L°° norm. This can be seen from Fig. 4, where the exponential decay 
of the L°° error in [0, oo[ between -D(A) and the cardinal expansion: 



D B (A) = J2 ^(A„)sinc 



n— — oo 



f(A-A n ) 
oc 



A n+ i — A n = Sc, 



(3.16) 



is displayed as a function of (5 C 1 on a Log-Lin plot Q. Switching back to the original variables, eq. ( 3.16 ) reads 

oo 

D B (M; 5 ' 3 -M T 5/3 ) := y D(M; 5 / 3 -M- 5 / 3 )-sinc ^(m t 5/3 -M~ 5 / 3 

On V 



where Q: 

M- 5 + { 3 -M- 5 ^ = S G . 
Given A4 S , the fitting factor obtained using ( |3.17 ) is given by: 



FF = maxD B (M-; 5/3 - M T 5/3 ) =: D B {M; 5/3 - M~ b/3 ) 



Mi 



It is convenient to let: 



so as to rewrite ( 3.1 9| ) as 



M: 5/3 = M- 5/3 + n6 c , »je[0,l[, 
M7 5/3 =M- 5/3 + V Jc, r?*G[0,l[, 

FF = D B (r)-r)*). 



(3.17) 

(3.18) 
(3.19) 

(3.20) 
(3.21) 

(3.22) 



The difference n — 77* turns out to depend on 7/ as shown in Fig. 5, and hence the fitting factor (3.22) depends on 77 
as shown in Fig. 6. The minimal-match condition (1.12) should again be enforced in the worst case(s), i.e., as seen 
from Fig.s 5, 6, for 77 = 0.5 and 77, = 77, yielding: 



71— — OQ 

















■ sine 




■i)-] 



r. 



(3.23) 



This condition is notably indipendent on q, and fixes the sample spacing Sc- 

In practice, as discussed in the previous section, it is convenient to use an economized cardinal expansion, viz: 

-5/3 



D B (M- 5/3 -M T 5/3 ):= D(M- 5/3 ~M- 5/3 )^ n (m^ /3 ~M^ 



n— — 00 



where (see eq.s (3.7) and ( 3.1C )): 



M^li 3 - M- 5 ' 3 = 5' c = x-'Sc 



(3.24) 



(3.25) 



and: 



sinh < 7rP(l — x 



v n (x) = sine ( — 



1 - 



PS'c 



1/2' 



PS' 



1/2 



(3.26) 



sinh[7rP(l-x^ 1 )] 
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As shown in the pr eviou s section, the function ( [3.26 ) is essentially contained in the interval: \x\ < P6' c , and hence 
the infinite sum in ( 3.24| ) is essentially restricted to: 



M. 



-5/3 



- P < n < 



M. 



-5/3 



P. 



Capitalizing on eq. ( 3.11 ) we shall enforce the condition: 

1 



sinh[7rP(l - x -1 ) 



i -r 

10 



(3.27) 



(3.28) 



to guarant ee th at the minimal-match con ditio n will not be affec ted w ithin the last significant fig ure o f I, w hen u sing 
using eq. (3.24) truncated according to (3.27) in place of eq. ( [3.17 ). Further, in view of eq.s ( j3"laj ) and fl3.15| ), we 
shall take: 



X = Xopt =■ 1 



26 c sinn" 1 



10 

i -r 



tt(m; 



5 /3_ 7 Q-5/3 



1/2 



(3.29) 



and: 



sinh 1 



P — Popt —'■ 



10 

i -r 



- Xopt) 



(3.30) 



so as to minimize the total number of correlators 



jvl min J vi -max ^ <^p 



(3.31) 



needed to evaluate (3.24) throughout the range [M m in, M. ma x\ of Mr- 



C. The Cardinal-Interpolated Reduced Correlator 



As a next step, we make the ansatz that an approximate representation of the reduced (newtonian) noncoherent 
correlator in terms of a (generalized) cardinal expansion also holds, viz. |46f| : 



C B := J2 C n y n (MT 5/3 - M- 5/3 ), 



(3.32) 



n— — oo 



where: 



max 

T 



f inf 



" A(f)f*(f) 

n(/) 



df 



(3.33) 



an d the infinite sum is truncated according to (3^). 

In (3.33) the templates T n are defined by eq.s Jl.9| ), (1.10) and (1.16), where the (scaled) chirp masses take the values 

M^ 5/3 = M~ 5 J 3 + kS' c , k = -P,-P + l,...,N c -l + P, (3.34) 

the interpolating functions ^/ n (x) are given by (|3.26| ), and the parameters 5c, X an d P are computed from the 
prescribed minimal match T as explained In Sect. III. A. 







IV. PLAIN VS. CARDINAL INTERPOLATED LATTICE 



In this section we shall compare the (newtonian) cardinal-interpolated (reduced, noncoherent) correlator to the 
plain-lattice of (reduced noncoherent) correlators in terms of computational cost and statis tical featu res (detection 
and estimation performance). The assumed noise PSD and spectral window are given by (1.21) and ( 1.22| ), respectively. 



A. Computational Burden 

The plain lattice template spacings 6l and total number of correlators Nl needed to cover the range (O.2M , 10M Q ) 
for some values of the minimal match T are compared in Table I to the corresponding quantities 8c and Nq of the 
cardinal interpolated correlator. 

It is seen that at any value of the minimal match T, the cardinal-interpolated representation requires some 30% 
less many templates than the plai n lat tice. 

On the other hand, evaluating ( |3.32| ) at any value Air ^ M-k is substantially cheaper than computing the corre- 
sponding (reduced , non coherent) correlator. 

Indeed, to use ( |3.32 ) one really needs to evaluate the interpolating functions ^ n {x) only at a finite number of 



(equispaced) values of A4 T 5 ^ 3 between the samples M. q 5 ^ 3 . The corresponding values of the interpolating functions 
can be computed once for all, and stored in a look-up table. As a result, only ~ 2P floating point operations are 
needed to compute (3.32) at each of the above values of 7W T 5 ^ 3 , with a typical P ~ 10 2 . 



B. Statistical Features 

The statistical properties of the cardinal-interpolated correlator have been compared to those of the plain-lattice 
via extensive Monte Carlo simulations. The number of different realizations used to derive the statistics was ~ 10 4 . 

Simulated data were sampled in time at twice the Nyquist rate. To limit running times, the minimal-match was 
set at r = 0.9, and the chirp-mass range was chosen in such a way that the longest observable waveform spanned 
2 15 time bins. In order to avoid circular-correlation artifacts, and to have equal statistics for all reduced correlators, 
all templates were zero-padded up to a total length of 2 16 bins. Gaussian uncorrelated noise samples were generated 
using a feedback-shift-register routine from the IMSL package, featuring an extremely large period |l7|], followed by 
a Box-Miiller transformation [Q. The noise samples were added to the whitened data in the spectral domain. 

In the following we shall denote the cardinal- interpolated and plain lattice test-statistics as M] : 



C£J X = ma^CA {Mt-Mk) (4.1) 

and: 



Mi 

k 



Ci L i=maxC fc . (4.2) 



respectively. We shall denote the estimated mass, i.e., the value of Mt which yields the maximum in (hi) or (4.2) 
as Ai es t, and let: 

M-* /3 =M- 5 / 3 +7 hst 8, rj est e[0,l[, (4.3) 

M-^^M-^z+vS, ?/e[0,l[, (4.4) 
Whe n ever needed a suffix/superfix L, C will be used to identify the plain-lattice and cardinal-interpolated cases in 



(4.3), (4.4) 



The CDFs of Cmax (dashed lines) and Cmax (full lines) in the absence of signal (SNR = 0), are compared in Fig. 
7-a, for template spacings corresponding to T — 0.9. The corresponding PDFs are displayed in Fig. 7-b. The observed 
difference falls within the 3tr uncertainty interval related to the finite number of realizations. 

The CDFs of Cmax (dashed lines) and Cmlx (full lines) in the presence of a signal with SNR — 6, 8, 10 are shown in 



Fig. 8-a, for the (worst) case where ry = 0.5 in (4.4). The corresponding PDFs are displayed in Fig. 8-b. Again, the 



observed differences fall within the 3cr uncertainty interval related to the finite number of realizations. Note that the 
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expected valu e alw ays exceeds the design value T ■ SNR, as might be expected as an effect of the supremum-taking 
operations in ( 1.15 ), (1.1) and (1.2). 

As r\ in (L4) changes between and 0.5, the PDFs of Cmax and Cmlx change in turn. The limiting PDFs 



corresponding to 77 = and 77 = 0.5 are shown in Fig.s 9 and 10, for the special case SNR = 8, for of C^L and C^} x . 

It can be concluded that the detection performance of the cardinal interpolated (reduced, noncoherent) correlator is 
essentially equivalent to that of the computationally more expensive plain lattice of (reduced, noncoherent) correlators. 
We turn now to a comparison of the pertinent estimation features. To this end we let: 



M~ 5 / 3 - M7 5/3 

6= ^ / " =V~Vesf (4.5) 



The PDFs of £ for the cardinal-interpolated correlator at SNR = 8 and T = 0.9 is shown in Fig. 11 for 77 = and 
rj = 0.5. 

The corresponding probabilities P(£) for the plain lattice of correlators are shown in Fig. 12. Note that for the 
plain lattice of correlators, rj est can take only values which are integer multiples of 8l- 

Both the cardinal-interpolated correlator and the plain lattice of correlators provide biased estimates. The bias 
E[£\ becomes for both essentially independent of the SNR at sufficiently high SNR levels (SNR ~ 8). For the 
cardinal-interpolated correlator the asymptotic large-SNR bias is shown in Fig. 13, which is a close akin of Fig. 5. 
For the plain lattice, it is displayed in Fig. 14. The cardinal-interpolated correlator is seen to exhibit a smaller bias. 

The std. deviations of the cardinal-interpolated and plain lattice estimators are nearly the same. For instance, 
for n = 0.5 (worst case), at SNR = 8 one has cr(£ L ) = 1.12 and cr(£ c ) = .92; at SNR = 10 a(£ L ) = 0.84 and 
<t(£c) = 0.63. 



C. Extension to PN Models 



The cardinal-interpolated approach can be extended in principle, to higher order PN models, provided the structure 
and q-BL properties of the (reduced) correlator are preserved. In the geometrical language of |f)|, this is equivalent 
to requiring that the chosen parameter space be (globally) flat and Euclidean [ 

models [fEiJ, and almost the case for the new 0— spin 2PN coordinates proposed in |3l). One should expect an even 
more substantial computational saving, in view of the higher dimension of the parameter space. 



50 1 . This is surely the case for 1PN 



V. CONCLUSIONS AND RECOMMENDATIONS 



Quasi-bandlimited function approximation theory can be used to build a (nearly) minimum redundant cardinal- 
interpolated representations of the noncoherent correlator for detecting gravitational wave chirps. An explicit expres- 
sion has been provided and tested, for the simplest case of newtonian waveforms. 

The number of correlators to be computed and interpolated in order to maintain the match above a given minimal 
value r has been shown to be substantially less then required by the std. (lattice) approach, and the computational 
gain goes up with T. On the other hand, evaluating the cardinal-interpolated representation at any value of A4t is 
substantially cheaper than computing the corresponding correlator. 

We suggest that cardinal-interpolated expansions could be used to improve the efficiency of hierarchical searches, 
at all hierarchical levels. 

Extension to PN templates should be straighforward, in principle, insofar as the structure and q-BL property of 
the correlators is preserved, and lead to an even increased computational gain. Work in this direction is in progress. 
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Fig. 1 - The function d(A,@). 




Fig. 2 - The function D(A), close-up. 
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Fig. 3 - Value of © which maximizes d(A,@) as a function of A (close-up). 




Fig. 4 - L error of cardinal expansion (3.16) versus sample density l/8c . 
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Fig. 5 - The difference r| - r|* as a function of r| - relevant to Eq.s (3.20)-(3.22). 




Fig. 6 - FF as a function of r|, relevant to eq.s (3.20), (3.22). 




Fig. 7a -Cardinal interpolated NCC (dashed) versus plain lattice of NCCs 
(continuous) at r=0.9. CDFs of detection statistic at SNR=0. 




Fig. 7b -Cardinal interpolated NCC (dashed) versus raw lattice of NCCs 
(continuous) at r=0.9. PDFs of detection statistic at SNR=0. 




Fig. 8a -Cardinal interpolated NCC (dashed) versus raw lattice of NCCs 
(continuous) at r=0.9, r|=0.5. CDFs of detection statistic at SNR=6,8,10. 



Lattice (cont.) 
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Fig. 8b -Cardinal interpolated NCC (dashed) versus raw lattice of NCCs 
(continuous) at r=0.9, r|=0.5. PDFs of detection statistic at SNR=6,8,10. 




Fig. 9 - Limiting PDFs of C { ^ max at SNR=8. 
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Table I - Plain versus cardinal interpolated lattice. 
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Fig. 10 - Limiting PDFs of C L; max at SNR=8. 




Fig. 1 1 - Cardinal-interpolated NCC at r=0.9. PDFs of estimated £ 
for r)=0 (dashed) and r|=0.5 (continuous), at SNR=8. 
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Fig. 12 - Plain lattice of NCCs at T=0.9. Probabilities of 
estimated £ for r\=0 (gray) and r|=0.5 (black), at SNR=8. 
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Fig. 13 - Cardinal interpolated NCC. Estimation bias (SNR > 8). 
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Fig. 14 - Plain lattice of NCCs. Estimation bias (SNR > 8). 



